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Abstract. The linear stability of MHD Taylor-Couette flow of infinite vertical extension is considered for various magnetic 
Prandtl numbers Pm. The calculations are performed for a wide gap container with 77 = 0.5 with an axial uniform magnetic 
field excluding counterrotating cylinders. For both hydrodynamically stable and unstable flows the magnetorotational instability 
produces characteristic minima of the Reynolds number for certain (low) magnetic field amplitudes and Pm > 0.01. For Pm 
< 1 there is a characteristic magnetic field amplitude beyond which the instability sets in in form of nonaxisymmetric spirals 
with the azimuthal number m — 1. Obviously, the magnetic field is able to excite nonaxisymmetric configurations despite of 
the tendency of differential rotation to favor axisymmetric magnetic fields which is known from the dynamo theory. If Pm is 
too big or too small, however, the axisymmetric mode with m=0 appears to be the most unstable one possessing the lowest 
Reynolds numbers - as it is also true for hydrodynamic Taylor-Couette flow or for very weak fields. That the most unstable 
mode for modest Pm proves to be nonaxisymmetric must be considered as a strong indication for the possibility of dynamo 
processes in connection with the magnetorotational instability. 
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1 . Introduction 

In order to discuss possible experimental realizations of the 
magnetorotational instability as the main transporter of angular 
momentum in all kinds of accretion disks there are several re- 
cent studies of Taylor-Couette flow for electro-conducting flu- 
ids between rotating cylinders under the influence of an uni- 
form axial magnetic field (Ji et al. 2001; Riidiger & Zhang 
2001; Willis & Barenghi 2002). The numbers describing the 
geometry of the container and the magnetic Prandtl number of 
the fluid have been considered as the free parameters. For a 
given magnetic field amplitude (the Hartmann number) the crit- 
ical angular velocity of the inner cylinder (the critical Reynolds 
number) is computed for the onset of an instabiUty of the rota- 
tion law between the cylinders. 

In Riidiger & Shalybkov (2002) the instability pattern is 
considered as axisymmetric. The main result for resting outer 
cylinder is that for high magnetic Prandtl number for weak 
magnetic field the excitation of the instability is easier than 
without magnetic field but for strong magnetic field the excita- 
tion of the instability is more complicated. The effect, however, 
disappears for small magnetic Prandtl number, i.e. for lower 
electric conductivity of the fluid as it may be realized in proto- 
planetary disks. 



On the other hand for rotating outer cylinder, when no insta- 
bility without magnetic field exists, the magnetic field always 
produces critical Reynolds numbers which, however, are run- 
ning with 1/Pm. For Pm of order 10^^ the critical Reynolds 
number is of order 10^ which is just the experimental limit. 

In the present paper the nonaxisymmetric perturbations are 
included into the consideration. This is of particular relevance 
for the question whether the Cowling theorem for dynamo ac- 
tion can be fulfilled, after which a dynamo can only work with 
nonaxisymmetric fields. We shall find that indeed for certain 
parameters - despite the smoothing action of the differential 
rotation - nonaxisymmetric modes can be excited easier than 
axisymmetric modes. This is in great contrast to earlier results 
of Taylor-Couette flow without magnetic fields where always 
the axisymmetric modes possess the lowest Reynolds numbers 
(Roberts 1965; DiPrima 1961)E|. 

Here, the dependence of a real Taylor-Couette flow on the 
magnetic Prandtl number and on the azimuthal 'quantum num- 
ber to' is investigated. The simple model of uniform density 
fluid contained between two vertically-infinite rotating cylin- 
ders is used with a constant magnetic field parallel to the ro- 
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' For counterrotating cylinders, however, the preference of nonax- 
isymmetric modes is already known, see Kriiger et al. 1966; Chen & 
Chang 1998 



2 



D.A. Shalybkov & G. Riidiger & M. Schultz: Nonaxisymmetric patterns of the linear MHD Taylor-Couette instability 





Olt 



Fig. 1. Cylinder geometry of the Taylor-Couette flow 



tation axis. The unperturbed state is a stationary circular flow 
with f] 



b/R^ 



(1) 



where a and b are two constants related to the angular velocities 
J7in and Q^ut with which the inner and the outer cylinders are 
rotating. If and i?out (^out > ^^in) the radii of the two 
cylinders then 
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After the Rayleigh stabihty criterion, d{R^QY /dR > 0, ro- 
tation laws are hydrodynamically stable for jl > ij^. Taylor- 
Couette flows with resting outer cylinders (jl = 0) are thus 
never stable. Here in order to isolate the MRI, we are also in- 
terested in flows with rotating outer cylinders so that the hydro- 
dynamical stability criterion, /i > 77^, is fulfilled. Our standard 
examples are formed with ij — 0.5, fi — and fi = 0.33, resp. 
The first example (fi = 0) is hydrodynamically unstable and 
the second one (fi = 0.33) is hydrodynamically stable. We are 
here only interested in flow patterns in containers with positive 

A- 



2. Basic equations 

The MHD equations which have to be solved are 
du 1 

— + (uV)u = --Vp + i'Au + J X B (4) 
ot p 

and 
dB 

^— = rot(M X B) + riAB (5) 

ot 

with divu = divB = 0. They are considered in cylindrical 
geometry with R, 0, and z as the cylindrical coordinates. A 
viscous electrically-conducting incompressible fluid between 
two rotating infinite cylinders in the presence of a uniform 
magnetic field parallel to the rotation axis leads to the ba- 
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sic solution Ub. 

const., and H — a + b/R^, with U as the flow and B as 
the magnetic field. We are interested in the stability of this so- 
lution. The perturbed state of the flow may be described by 



Us 



Bq + B'^, p', with p' as the 
pressure perturbation. 

Here only the linear stability problem is considered. By an- 
alyzing the disturbances into normal modes the solutions of the 
linearized magnetohydrodynamical equations are of the form 

u'^ = M^(i?)e'(™'*+'=^-'^*), B'^ = B0(i?)e'('"'^+'=^-"*\ 

u'^ = u^(i?)e'(™'^+'=^-'^'), B'^ = B,{R)e'^""l'+'"'-'^'\ (6) 

Only marginal stability will be considered where the imaginary 
part of uj vanishes. Let d = i?out — Rin be the gap between the 
cylinders. We use 

(i?ind)'/' (7) 



H 



as unit of length, the rj/H as unit of perturbed velocity and Bq 
as unit of perturbed magnetic field with the magnetic Prandtl 
number 



Pm = 



V 



(8) 



1/ is the kinematic viscosity, rj is the magnetic diffusivity. Note 
H^^ as the unit of wave numbers and ly/H^ as the unit of 
frequencies. After elimination of both pressure fluctuations and 
the fluctuations of the vertical magnetic field, B'^, the equations 
are 

dR R ' R 



-fT + -f^u^ + \kUz — 0, 
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The Reynolds number Re and the Hartmann number Ha are 
defined as 



Re = 



Ha: 



B^H 



(14) 



For given Hartmann number and Prandtl number in the present 
paper we shall derive in a linear theory the critical Reynolds 
number of the rotation of the inner cylinder for various mode 
numbers m. 

3. Boundary conditions, numerics 

An appropriate set of ten boundary conditions is needed to 
solve the system (P)..([l3|). Always no-slip conditions for the 
velocity on the walls are used, i.e. ur = — duR/dR — 0. 
The boundary conditions depend on the electrical properties of 
the walls. The transverse currents and the perpendicular com- 
ponent of the magnetic field vanish on conducting walls hence 
dB^ I dR + i?0 /R — Br — Q. These boundary conditions hold 
both for R = i?;^ and for R = i?out- 

The homogeneous set of equations (p|)...(|l3|) with the 
boundary conditions for conducting walls determine the eigen- 
value problem of the form m, Re, Ha, 7?.((jj)) = for 
given Pm. The real part of uj, TZiuj), describes a drift of the pat- 
tern along the azimuth which only exists for nonaxisymmetric 
flows. £ is a complex quantity, both its real part and its imagi- 
nary part must vanish for the critical Reynolds number (Fig. ||). 
The latter is minimized by choice of the wave number k. TI{lo) 
is the second quantity which is fixed by the eigenequation. 

The system is approximated by finite differences with typ- 
ically 81 gridpoints. The resulting determinant, C, takes the 
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Fig. 2. The zero-lines of the real (solid) and the imaginary 
(dashed) part of the system determinant. At the crossing point 
both parts of the determinant simultaneously vanish fixing the 
Reynolds number and the drift frequency TI{lo). Pm and kz are 
prescribed 



value zero if and only if the values Re are the eigenvalues. For 
a fixed Hartmann number, a fixed Prandtl number and a given 
vertical wave number we find the eigenvalues of the equation 
system. They are always minimal for a certain wave number 
which by itself defines the marginally unstable mode. The cor- 
responding eigenvalue is the desired Reynolds number. 

4. Resuits for conducting waiis 

Only a container is considered in the present paper with one 
and the same gap geometry, i.e. = 0.5. Then the flow between 
the cylinders is hydrodynamically unstable between /i = and 
/t = 0.25. We shall work with both the hydrodynamically un- 
stable container with /i = and with the hydrodynamically 
stable container with /i = 0.33. 

4. 1. Resting outer cylinder (steep rotation law) 

We start with the results for containers with resting outer cylin- 
ders (Fig. 1^. Provided a critical rotation rate of the inner cylin- 
der is exceeded they are hydrodynamically unstable. Of course, 
for Ha = m = the known critical Reynolds number Re = 68 
is reproduced. For m > the critical Reynolds numbers ex- 
ceed the value for m — 0. Without magnetic field the instabil- 
ity yields rolls. The critical Reynolds number for m = 1 is 73 
and for ?7i = 2 it is 101 (Roberts 1965). 

With magnetic fields (Ha >0) the magnetic Prandtl number 
comes into the game. Results for Pm = 10, 1, 0.1 and 0.01 are 
presented in Fig. ||. For Pm > 1 the electrical conductivity is so 
high that the magnetorotational instability (Balbus & Hawley 
1 99 1 ; Brandenburg et al. 1 995 ; Ziegler & Rudiger 2000) for Ha 
~ 5 produces a characteristic minimum of the critical Reynolds 
number but for stronger magnetic fields the suppressing action 
of the magnetic field starts to dominate. In contrast to the ex- 
pectations, however, for the magnetic Prandtl numbers which 
are not too high and not too low the mode with m = 1 be- 
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comes more and more dominant. This is a new and interesting 
result: The linear instability of the Taylor-Couette flow without 
magnetic field is formed by axisymmetric rolls but the mag- 
netic field favors the excitation of bisymmetric spirals. For Ha 
> 10 ... 20 the instability sets in in form of a drifting pattern 
with maximum and minimum separated by 180°. However, as 
can be seen in Fig. ^(last plot) for small magnetic Prandtl num- 
ber (here Pm=0.01) again the axisymmetric pattern with m = 
again starts to dominate with the lowest critical Reynolds num- 
ber. 



The modes with ni = 2, which we have also computed, do 
never possess the lowest Reynolds numbers, they are not im- 
portant for the discussion of the pattern of the instability. What 
we have found is that in contrast to the hydrodynamic case 
(Pm=0) there are experimental combinations where the non- 
axisymmetric mode with m = 1 has a lower Reynolds number 
than the axisymmetric mode with to = 0. This is one of the 
most surprising structure-forming consequences of the inclu- 
sion of magnetic fields to the Taylor-Couette flow experiment 
found first in astrophysical simulations. 




Fig. 3. Resting outer cylinder: Stability lines for axisymmetric 
(to = 0, solid lines) and nonaxisymmetric instability modes 
(to = 1 (dashed lines), m ~ 2 (dashed-dotted lines). Results 
are given for Pm — 10, 1, 0.1 and 0.01. Note that for Pm — 1 
and for Pm = 0.1 for certain magnetic fields the nonaxisym- 
metric modes with to = 1 possess the lowest Reynolds num- 
bers 



4.2. Rotating outer cylinder (flat rotation law) 

If the outer cylinder rotates with an angular velocity fi > ff' 
than the linear instability without magnetic field disappears and 
the critical Reynolds number for Ha = moves to infinity. 
However, for finite Hartmann number (again of order 10) the 
instability survives practically for the same Reynolds numbers. 
The consequence is the occurrence of typical minima in the 
stability diagram (Fig. ^for/i = 0.33). 

The minima also occur for the nonaxisymmetric solutions 
with TO = 1. For very high electrical conductivity (Pm=10) 
there seems to be no intersection between both the bifurcation 
profiles. The ring-like structure with to = always possesses 
the lowest critical Reynolds number. 

This is not true, however, for smaller magnetic Prandtl 
numbers, i.e. for lower electrical conductivity. For Pm < 1 we 
always find intersections between the lines for to = and 
m — 1. Again there is a critical Hartmann number at which 
the ring geometry (to — 0) of the excited flow and field pattern 
changes to a nonaxisymmetric geometry with to — 1.^ 

Hence, also in experiments with rotating outer cylinder the 
magnetic field is able to produce nonaxisymmetric structures. 
After the Cowling theorem which requires the existence of non- 
axisymmetric magnetic modes for the existence of a dynamo, a 
selfexcited dynamo might thus exist, but only for certain mag- 
netic Prandtl numbers, i.e. for Pm < 1. The magnetic Prandtl 
number for experiments with liquid metals like sodium or gal- 
lium with Pm of order 10^'^ are still smaller than the con- 
sidered values. 

5. Wave number and drift frequencies 

The wave numbers have been discussed in detail for the ax- 
isymmetric modes in a foregoing paper (Riidiger & Shalybkov 
2002). Generally, the cells for the nonaxisymmetric modes be- 
come more and more elongated in the vertical direction. Here 
we only add remarks about the drift velocity TZ{u!) which al- 
ways proved to be positive, i.e. the pattern drifts are in direction 
of the rotation (eastward). It is 

• 7^(a;)^?in 

(P^ 5 , (15) 

TOKe 

so that for to = 1 the drift period in units of the rotation period 
runs as RdTZ{uj). A typical value for this ratio is 2. In Fig. ^ 

^ For differentially rotating spheres we find a similar behavior 
but only for stress-free boundary conditions (Kitchatinov & Riidiger 
1997) 
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Fig. 4. The same as in Fig. |^ but for rotating outer cylinder 
(fi — 0.33), there is no hydrodynamical instability 

the expression (|l5|) is given normalized with the rotation rate 
J?in for the solutions with m — 1 and for the two low magnetic 
Prandtl numbers that we have considered. 



Pm=0.1 Pm=0.01 
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Fig. 5. The drift frequencies normalized with the rotation 
rate J?in for the solutions with m = 1 for resting outer cylin- 
ders. The drift of the spirals is always positive, i.e. in the direc- 
tion of the rotation with about 50 % of the rotation of the inner 
cylinder 

On the other hand, after the Cowling theorem of dynamo 
theory a magnetic field can only be maintained if it is nonax- 
isymmetric. Considering a number of typical magnetic Prandtl 
numbers we find that for our container with conducting cylin- 
ders the dominance of the nonaxisymmetric modes only oc- 
curs for not too high and not too low magnetic Prandtl num- 
ber Obviously, the dissipation processes are more important 
for nonaxisymmetric modes rather than axisymmetric modes. 
Hence the dissipation allows nonaxisymmetric modes only to 
be preferred if both the dissipation values have nearly the same 
order of magnitude. 
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6. Conclusions 

We have shown that a Taylor-Couette flow which is stable in the 
hydrodynamic regime (p, > ff) is destabilized by a weak ax- 
ial magnetic field. Below a critical Hartmann number of order 
10... 100 the instability sets in in form of axisymmetric rolls 
while above this value the instability forms nonaxisymmetric 
field and flow modes. This phenomenon exists despite of the 
observation (e.g. in dynamo theory) that differential rotation is 
known as suppressing the formation of nonaxisymmetric mag- 
netic fields. 



